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Abstract. The theory of gauge-invariant non-spherical metric perturbations of Schwarzschild 
black hole spacetimes is now well established. Yet, as different notations and conventions have 
been used throughout the years, the literature on the subject is often confusing and sometimes 
confused. The purpose of this paper is to review and collect the relevant expressions related 
to the Regge- Wheeler and Zerilli equations for the odd and even-parity perturbations of a 
Schwarzschild spacetime. Special attention is paid to the form they assume in the presence of 
matter-sources and, for the two most popular conventions in the literature, to the asymptotic 
expressions and gravitational-wave amplitudes. Besides pointing out some inconsistencies in 
the literature, the expressions collected here could serve as a quick reference for the calcula- 
tion of the perturbations of Schwarzschild black hole spacetimes driven by generic sources and 
for those approaches in which gravitational waves are extracted from numerically generated 
spacetimes. 



PACS numbers: 04.30.Db, 04.40.Dg, 04.70.Bw, 

1. Introduction 

Black-hole perturbation theory (see for example the monograph by Chandrasekhar [ 1 1, the 
book by Frolov and Novikov [2 1 or the recent review by Kokkotas and Schmidt 1 3 1) has been 
fundamental not only for understanding the stability and oscillations properties of black hole 
spacetimes |4), but also as an essential tool for clarifying the dynamics that accompanies 
the process of black hole formation as a result of gravitational collapse |6). As one 
example among the many possible, the use of perturbation theory has led to the discovery that 
Schwarzschild black holes are characterised by decaying modes of oscillation that depend on 
the black hole mass only, i.e. the black hole quasi-normal modes 171 151 l9l flOl . Similarly, 
black-hole perturbation theory and the identification of a power-law decay in the late-time 
dynamics of generic black-hole perturbations has led to important theorems, such as the 'no 
hair' theorem, underlining the basic black-hole property of removing all perturbations so that 
'all that can be radiated away is radiated away' p||6| ll II . 

The foundations of non-spherical metric perturbations of Schwarzschild black holes 
date back to the work in 1957 of Regge and Wheeler |4|, where the linear stability of the 
Schwarzschild singularity was first addressed. A number of investigations, both gauge- 
invariant and not, then followed in the 1970s, when many different approaches were proposed 
and some of the most important results about the physics of perturbed spherical and rotating 
black holes established 0|6]|3|![To|[I2[T3[l![l3[T! Building on these studies, 

which defined most of the mathematical apparatus behind generic perturbations of black 
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holes, a number of applications have been performed to study, for instance, the evolutions 
of perturbations on a collapsing background spacetime fT9l [20l f2Tl f22l [231 f24l [251 [26l. 
Furthermore, the gauge-invariant and coordinate independent formalism for perturbations 
of spherically symmetric spectimes developed in the 70's by Gerlach and Sengupta in 
1191 1201 |2"T1 . has been recently extended to higher-dimensional spacetimes with a maximally 
symmetric subspace in 1271 1281 l29l 1301 . for the study of perturbations in brane-world models. 

Also nowadays, when numerical relativity calculations allow to evolve the Einstein 
equations in the absence of symmetries and in fully nonlinear regimes, black hole perturbative 
techniques represent important tools. Schwarzschild perturbation theory, for instance, has 
been useful in studying the late-time behaviour of the coalescence of compact binaries in 
a numerical simulation after the apparent horizon has formed H3 1 1 1321 l33l . In addition, 
methods have been developed that match a fully numerical and three-dimensional Cauchy 
solution of Einstein's equations on spacelike hypersurfaces with a perturbative solution in a 
region where the components of three-metric (or of the extrinsic curvature) can be treated as 
linear perturbations of a Schwarzschild black hole (this is usually referred to as the 'Cauchy- 
Perturbative matching') 1341 1351 l36l 1371 1381 . This method, in turn, allows to "extract" the 
gravitational waves generated by the simulation, evolve them out to the wave-zone where 
they assume their asymptotic form, and ultimately provide outer boundary conditions for the 
numerical evolution. 

Overall, therefore, black hole perturbation theory represents a very powerful tool with 
multiple applications and developments. The purpose of this review is to review the 
theory of metric perturbations of Schwarzschild black holes, especially in its gauge-invariant 
formulations. While a systematic derivation of the most relevant expressions is provided, 
special care is paid to 'filter' those technical details that may obscure the important results 
and provide the reader with a set of expressions that can be readily used for the calculation 
of the odd and even-parity perturbations of a Schwarzschild spacetime in the presence of 
generic matter-sources. Also, an effort is made to 'steer' the reader through the numerous 
conventions and notations that have accompanied the development of the formalism over 
the years, pointing out the misprints and errors that can be found in the literature. Finally, 
although most of what presented here is not new, some expressions that have not been 
discussed before in the literature are also introduced. 

It should be noted that our discussion will be restricted exclusively to black-hole 
spacetimes, although a Schwarzschild metric can also be used to describe the vacuum exterior 
of a spherical relativistic star. The reason behind this choice is that boundary conditions 
different from the ones used for black holes need to be imposed at the stellar surface. 
Furthermore, in the case of a relativistic star the Einstein equations need to be suitably 
coupled with the equations of relativistic hydrodynamics. A detailed discussion of how 
to do this in practice is beyond the purpose of this review, which is instead focused on a 
vacuum background spacetime. However, a recent review on the theory of relativistic stellar 
perturbations can be found in 1 3 1. 

The organisation of the paper is as follows. In section|2]we review the basic properties 
and advantages of a gauge-invariant description of the perturbations of a curved spacetime, 
while in section|3]we review the multipolar expansion of the perturbations of a Schwarzschild 
spacetime. Sections 0] and [5] are instead devoted to the derivation of the general equations 
governing the evolution of odd and even-parity perturbations, respectively, whose asymptotic 
expressions for the gravitational-wave amplitudes and losses are presented in section [6] 
Finally, section Q contains our conclusions and a number of Appendices follow to provide 
some useful relations and expressions for the angular parts of the perturbations. 

We use a spacelike signature (—,+,+,+) with Greek indices running from to 3 and 
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indicate partial and covariant derivatives in the /i-th direction with the symbols and V M , 
respectively. Also, symmetrized and antisymmetrized indices in a tensor will be indicated 
with round and square brackets, respectively. 

2. Gauge-invariant metric perturbations 

To highlight the importance of gauge-invariant quantities in describing the metric 
perturbations it is useful to recall that even if the coordinate system of the background 
spacetime has been fixed, the coordinate freedom of general relativity introduces a problem 
when linear perturbations are added. In particular, it is not possible to distinguish an 
infinitesimal 'physical' perturbation from one produced as a result of an infinitesimal 
coordinate transformation (or gauge-transformation). This difficulty, however, can be 
removed either by explicitly fixing a gauge (as done in the pioneering works on the 
subject 0] 15] |6] ID [T2l ElX or by introducing linearly gauge-invariant perturbations (as 
initially suggested by Moncrief 1 14 1 and subsequently adopted in several applications HI 5111 61 

ESlEUESll^l). 

More specifically, given a tensor field X in a background metric g and SX its 
infinitesimal perturbation, an infinitesimal coordinate transformation x^ — > x M = x^ + £ M 
with £ p <C 1 will yield a new tensor field 

5X -> SX' = SX + C^X , (1) 

where is the Lie derivative along £ in the metric g. We will then consider SX to be 
gauge-invariant if and only if C^X = 0. 

Stated differently, the possibility of building gauge-invariant metric perturbations relies 
on the existence of symmetries of the background metric. In the case of a general spherically 
symmetric background spacetime (i.e. one allowing for a time dependence) and which has 
been decomposed in multipoles (see Sect. [3}, the construction of gauge-invariant quantities 
is possible for multipoles of order £ > 2 only 1191 1201 1391 1401 . In practice, the advantage 
in the use of gauge-invariant quantities is that they are naturally related to scalar observables 
and, for what is relevant here, to the the energy and momentum of gravitational waves. At the 
same time, this choice guarantees that possible gauge-dependent contributions are excluded 
by construction. 

Of course, this procedure is possible if and only if the background metric has the 
proper symmetries under infinitesimal coordinates transformation and a gauge-invariant 
formulation of the Einstein equations for the perturbations of a general spacetime is not 
possible. Nevertheless, since any asymptotically flat spacetime can in general be matched 
to a Schwarzschild one at sufficiently large distances, a gauge-invariant formulation can be 
an effective tool to extract physically information about the gravitational waves generated in 
a numerically evolved, asymptotically flat spacetime 1341 1351 l36ll37M38l . 

3. Multipolar expansion of metric perturbations 

Given a spherically symmetric Schwarzschild solution with metric g and line element 

ds 2 = g^dx^dx" = -e 2a dt 2 + e 2b dr 2 + r 2 (d9 2 + sin 2 6d<j) 2 ) , (2) 

where e 2a = e~ 2b = (1 — 2M/r), we generically consider small non-spherical perturbations 
h^v such that the new perturbed metric is 

9i±v — 9 i±i> ~l~ h^iv j (3) 
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where <C 1- It should be noted that although we have here chosen to employ 

Schwarzschild coordinates to facilitate the comparison with much of the previous literature, 
this is not the only possible choice. Indeed, it is possible to formulate the perturbations 
equations independently of the choice of coordinates as discussed in |49| and in [51 1. 

Because the background manifold M. is spatially spherically symmetric, it can be 
written as the product M. = M 2 x S 2 , where M 2 is a Lorentzian 2-dimensional manifold 
of coordinates (t,r) and S 2 is the 2-sphere of unit radius and coordinates (9,(f>). As a 
result of this decomposition, the perturbations can be split ab initio in a part confined to 
M 2 and in a part confined on the 2-sphere S 2 of metric 7. Exploiting this, we can expand 
the metric perturbations h in multipoles referred to as odd or even-parity according to their 
transformation properties under parity. In particular, are odd (or axial) multipoles those that 
transform as (— 1) +1 , under a parity transformation (9, </>) — > (tt — 9, tt + </>), while are even 
(or polar) those multipoles that transform as (— l) e . As a result, the metric perturbations can 
be written as 

■£ra\(°) , /ilra\( e ) 



£,m 



[IIS } 



(4) 



where 



5>E 
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Introducing now a notation inspired by that of Gerlach and Sengupta 1191 1201 1211 and 
recently revived by Gundlach and Martin-Garcia 11391 140| 1411 . we use upper-case indices 
A, B . . . — 0, 1 to label the coordinates of M 2 and lower-case indices c,d. . . = 2, 3 to 
label the coordinates of S 2 . 

Using this notation, the scalar spherical harmonics Y lm are then simply defined as 



7 



cd 



irri 



-AY 



I'm 



(5) 

where V c indicates the covariant derivative with respect to the metric 7 = diag(l, sin 2 9) of 
S 2 and where A = £(£ + 1). It is now convenient to express the odd and even-parity metric 
functions in in terms of tensor spherical harmonics. To do this we introduce the axial 
vector S^ m defined as 

S l c m = e cdl de W e Y em , (6) 

where e c d is the volume form on S 2 as defined by the condition e C( je ce = 7 d e and such that 
V c e a b = 0. In this way, each odd-parity metric function in Q can be written as 



(o) 



( 



(7) 



where , h are functions of (t, r) only and where we have omitted the indices £, m for 
clarity. 

Proceeding in a similar manner, each even-parity metric function can be decomposed in 
tensor spherical harmonics as 



/ e 2a H Q Y lm Hi 



tm 



H{Y im e 2b H 2 Y lm 



V h^V c Y em 



h ( ^y c Y em 



(8) 



r 2 (KY lm lcd + GV d V c Y lm ) J 
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where H , Hi, H 2 , K, G (with the indices £, m omitted for clarity) are the 

coefficients of the even-parity expansion, are also functions of (t, r) only. 

Note that we have here used the Regge-Wheeler set of tensor harmonics to decompose 
the even-parity part of the metric in multipoles [4|. Despite this being a popular choice in 
the literature, it is not the most convenient one since the tensor harmonics in this set are 
not linearly independent. An orthonormal set is instead given by the Zerilli-Mathews tensor 
harmonics 1421 1431 and the transformation from one basis to the other is given by defining the 
tensor Z^J 1 confined on the 2-sphere S 2 



/y£rn WW \^im 



A 



:lc d Y 



Cm 



(9) 



and then replacing in equation (|8j the second covariant derivative of the spherical harmonics 
V C V dX lm with Z^J 1 . This transformation has to be taken into account, for instance, when 
developing gauge-invariant procedures for extracting the gravitational-wave content from 
numerically generated spacetimes which are "almost" Schwarzschild spacetime 1441 1451 l46l 
|47]|48). 

Besides vacuum tensor perturbations, the background Schwarzschild spacetime can be 
modified if non-vacuum tensor perturbations are present and have a mass-energy much smaller 
than that of the black hole but sufficiently large to perturb it. In this case, the generic stress- 
energy tensor t^ u describing the matter-sources can be similarly decomposed in odd and even- 
parity parts 



v=E[(a (0) 



(o (c) 



(10) 
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that are naturally gauge-invariant since the background is the vacuum Schwarzschild 
spacetime and are given explicitely by 

i gin 

«™) vo) =( I <■» 



j lm at 



for the odd-parity part and by 

(a (c) - 



1 AB 1 



\ T| m v c y 



em 



T^ % V c Y tm 



(12) 



r J 3 r jcd + J- 2 ^cd I 



for the even-parity one. Note that we have now used the Zerilli-Matthews set of harmonics for 



the expansion, that the ten coefficients L^ 71 , L , Tjjg, T^ 7 



T2" 1 , T| m are gauge-invariant 



and that explicit expressions for them will be presented in the following Sections. 

With the perturbations decomposed in a convenient form, we can now consider the 
Einstein field equations that, in the static vacuum background, take the simple form 



Rnu — 



(13) 



where R is the Ricci tensor built from the background metric g. At first order in the 
perturbations, the linearity can be exploited to break up the Einstein equations as 



1? 



(14) 
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where R is now the Ricci tensor built from the metric perturbations h. Using equations dl 31 . 
the field equations then reduce to 

Flfiv - x9fiuR = Sirt^v , (15) 
and where the components of the Ricci tensor are explicitely given by 

R»u = - v„rjg . (16) 

Here, T are the Christoffel symbols relative to the perturbed metric h, i.e. 



and should be distinguished from the corresponding Christoffel symbols T relative to the 
background metric g. 

Note that while a generic perturbation will be a mixture of odd and even-parity 
contributions, we will exploit the linearity of the approach to handle them separately and 
simplify the treatment. In the following two Sections we will discuss the form the Einstein 
equations (I15> assume in response to odd and even-parity perturbations over a Schwarzschild 
background. In particular, we will show how the three odd-parity coefficients of the expansion 
in harmonics of the metric and the seven even-parity ones can be combined to give two gauge- 
invariant master equations, named after Regge and Wheeler 1 4 1 and Zerilli 1 42 1 , each of which 
is a wave-like equation in a scattering potential!. 

Although our attention is here focussed on the radiative degrees of freedom of the 
perturbations (i.e. those with t > 2) because of their obvious application to the modelling of 
sources of gravitational waves, a comment should be made also on lower-order multipoles. 
In particular, it is worth remarking that the monopole component of the metric for a vacuum 
perturbation (i.e. with I = Q) is only of even-parity type and represents a variation in the 
mass-parameter of the Schwarzschild solution. On the other hand, the dipole component of 
the even-parity metric for a vacuum perturbation (i.e. with I = 1) is of pure-gauge type 
and it can be removed by means of a suitable gauge transformation 111 31 . Note that this is 
not the case for a dipolar odd-parity metric perturbation, which can instead be associated to 
the introduction of angular momentum onto the background metric. A detailed discussion of 
these low multipoles can be found in Appendix G of 1 13 1 and in |49 1. 



4. Gauge-invariant odd-parity perturbations 

Before discussing the derivation of the odd-parity equation a choice should be made for the 
odd-parity master function. Unfortunately, this choice has not been unique over the years 
and different authors have made different choices, making comparisons between different 
approaches less straightforward. Here, we will make a choice which highlights the relation 
with the gravitational-wave amplitude measured by a distant observer and, in particular, we 
construct the gauge-invariant combination of multipoles [ 19 39"! 1501 

Vat 

k A = h A - V A h + 2h—^- , (18) 

r 

t. These results were originally obtained by Regge and Wheeler |4 and by Zerilli 1 1 21 1131 in a specific gauge 
(i.e. the Regge-Wheeler gauge). Subsequently, the work of Moncrief showed how to reformulate the problem in a 
gauge-invariant form by deriving the equations from a suitable variational principle 1 14 . Some details of Moncrief 's 
approach to odd-parity perturbations are briefly described in |Appendix C| 
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where, we recall, represents the covariant derivative with respect to the connection of the 
submanifold M 2 . In this way, the function 



$<°>(t, r) 



r 3 e AB 



(o) 



r 2 d r 



(19) 



where eAB is the antisymmetric volume form on M 2 (e 01 = — e 10 = — 1, see| Appendix Q , is 
gauge-invariant and will be our choice for the Regge- Wheeler master function J19II39II40II50I . 

A slight variation of the master function dl9t has been introduced by Cunningham, Price 
and Moncrief 1 15 1 in terms of the function if) = A$'°' and this has been used so extensively 
in the literature 1231 l24l f'26 1 that it is now commonly referred to as the Cunningham-Price- 
Moncrief (CPM) convention. We partly follow this suggestion and introduce a new master 
function for the odd-parity perturbations defined as 

\$r(°) = __L_ <j>(°) . (20) 

With the choice ( I20K the Einstein field equations dl 51 with odd-parity perturbations lead 
to the inhomogeneous 'Regge- Wheeler' equation 

<9 2 * (o) - d 2 r ^ {o) + V e (o) ^ = S (o) , (21) 

where 



2Mln 



2M 



- 1 



is the 'tortoise coordinate' II II and v} is the odd-parity potential, defined as 

A 6M\ 



(22) 



(23) 



The right-hand-side of equation ( 12 1 1 > represents the generic odd-parity "matter-source" and is 
given by 



S (o) 



167T?' 



„2aAB 



167rr 



1 



2M 
r 



d t L\ m -d r ,L 



A-2 A-2 
with the components of the odd-parity matter-source vector defined as 



t £m 



sinf 



(i m t A2 Y e * n + t A3 d e Y t * r 



im) 



A = 0,1 



(24) 



(25) 



and where dil = sin 6d(f)d6 is the surface element on the 2-sphere S 2 . 

As a result of this construction, given a generic odd-parity matter-source whose 
dynamics is known (e.g. through the solution of the equations of motion if the source is point- 
like, or those of relativistic hydrodynamics if the source is extended), equation (12 It describes 
the evolution of the perturbations induced on the black-hole spacetime by the dynamics of the 
matter-source. 

It should be noted that another choice for the gauge-invariant odd-parity master variable 
is possible and indeed was originally proposed by Moncrief [14|. This function, which 
hereafter we will refer to as the odd-parity Moncrief function, is defined as 



)(o) 



Vat oAb 
9 k B 



1 



,(°) 



-Or 



(26) 



where the first expression is coordinate independent |49|, while the second one is specialized 
to Schwarzschild coordinates with h 2 = —2h [14|. In the Regge- Wheeler gauge h 2 = 
h = 0, and the definition (I26> coincides with the variable used by Regge and Wheeler 1 4 1 . 
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Historically, the choice of J26I as master variable has been the most common in the literature 
to describe odd-parity perturbations of a Schwarzschild spacetime and we will refer to it 
as 'Regge-Wheeler' (RW) convention. It should be noted that while J26> is a solution of 
the Regge-Wheeler equation, the corresponding source term differs from expression J24l >. A 
general expression of the source in the RW convention can be found in |49 52 1 together with 
its specification for a point-particle (see also I53II54I 55)). 

The two master functions Q(°' and are intimately related in at least two different 
ways. Firstly, through the variational formalism employed by Moncrief in 1 14| and a brief 
discussion of these relations is presented in |Appendix C| Secondly, through an explicit 
expression that reads |49 1 

^(o) = _Q(o) + i^_^_ L ^. (2?) 

Equation J27> highlights an important difference between the two master functions which 
is not just a dimensional one (i.e. has the dimensions of a length, while Q(°> 

is dimensionless) and that will have consequences on the asymptotic expressions for the 
gravitational waveforms when these are expressed in one or in the other convention. A detailed 
discussion of this will be made in sections. 16. Ill6.3landl6.4l 



5. Gauge-invariant even-parity perturbations 

Also in the case of even-parity perturbations, an evolution equation similar to the Regge- 
Wheeler one i2l\ can be found. In particular, following Moncrief 1 14 1, we define the functions 

ki = K + (rd r G - h[ c) ^j , (28) 

k 2 = I [e 2b H 2 - e% (re b K)} , (29) 

where both n\ and k 2 are gauge-invariant functions, as well as the following linear 
combination 

4r 

q 1 = rA«i + -rrn 2 . (30) 



Strictly related to expression (I30i is the gauge-invariant function most frequently used in the 
literature ll231l4T1l52ll5^l57ll5^l59ll^0ll62 IS 



*W = V M (31) 

" A [r (A - 2) + 6M] V ; 

which is also the solution of the inhomogeneous even-parity master equation or 'Zerilli' 
equation 

<9 t 2 * (c) - 9 r 2 t * (c) + ^ (c) * (c) = , (32) 
and, again, is a wave-like equation in the scattering Zerilli potential 1121 



2M \ A(A - 2) 2 r 3 + 6(A - 2) 2 Mr 2 + 36(A - 2)M 2 r + 72M 3 
~ ~) r 3 [(A-2)r + 6M] 2 ' < " 1 ' 1 ' 

§ Note that expression 1311 con'ects the sign in the corresponding definition of the even-parity function used in 1 63 
(cf equation(7) in 1631 1. 
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The even-parity matter-source has a rather extended expression given by 



5 (e) = _ 



A( 6r 3 - WMr 2 ) - r 3 A 2 - 8r 3 + 68Mr 2 - 108Af 2 r 



A [(A- 



2)r + 6M] 
1 



rAl> 



e 

4Ar 



2Afr + r 2 (A-4) 



-n t r 

1 11 



(A - 2)r + 6M 
2r 3 



2r 3 9 r T n ^™ 



„4b 



Q rplr, 
U T» ± 11 



rpirm | ' 



2A 1 



A 



T 




" J 00 



(34) 



and we refer the interested reader to [52 60 , 61 1 where the steps necessary for the derivation 
of equation (I34> are given in detail and some application to accretion problems are discussed. 
Here, we simply recall that the expressions of the even-parity vector and tensor spherical- 
harmonics for the matter-source needed in J34i can be obtained from the orthogonality 
properties of the harmonics and are 



T>£rn 
A 



A 



(19. 



t A 2de{Y e * m ) - t A3 



vrnYl 



dm 



T. 



rptm 

J 3 — 



A(A 

1 

2T 2 " 



^22" 



w, 



t III 



dfl <22 + t 3 3^-o- 

sm i 



+ t 



sin 2 9 
2X 



23" 



irn 



+ t. 



X>, 



AY, 



I'm 



A = 0,1, 
m 2 Y* 



(35) 



cot OdeYZ 



sm 



Y 



t ab — J dCl t AB Y g * m , 
where the angular functions W 



£m ) 



and AT* 



W 



X 1 



diY 



sm f 



■sin0 V fl 5|" 



- cot 6d e Y em 



A, .0 = 0,1 
<p) are defined as |4| 

:dlY lm , 



sm 



= 2{d 2 H Y lm -cotBd^Y 1 



(36) 
(37) 

(38) 



(39) 
(40) 



and where the asterisk stands for complex conjugation. A few comments are now needed 
about the different notations in which the even-parity functions can be found in the literature. 
A particularly common choice is that proposed by Moncrief in 1 14 1 for the the even-parity 
gauge-invariant master function which is related to Zerilli function ( 13 U simply as 
Q(<0 = Af( c ). Other authors (e.g. (52] |56)) use instead a master function defined as 
Z = 2^( c ) and a factor of 2 needs to be introduced or removed when converting between 
the two notations. Yet another even-parity function can be introduced in terms of two new 
gauge-invariant metric functions k and \ where I 40ll41l 
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= K + 
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k = 


H 2 - 
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72h 



rd r G--h\ 
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(e) 



2M, 



(e) 



±- b dr(r 2 d r G) 



Md r G 



and such that 



k 2 



„2b 



X — rd r k + 



M 



2b, 



(41) 
(42) 

(43) 
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In this case, the Zerilli function ( 13 1 1 can be equivalently defined as |41 1 



\]/( e ) 



2r 2 



M 



,2li 



k — rd r k 



(44) 



A [(A - 2)r + 6M] e 2b 

Finally, we recall that the homogeneous odd and even-parity master equations (12 It and 
d32l > can be transformed into each other by means of differential operations 1 1 1 and they are 
connected to the master equation that Bardeen and Press |65 1 have derived via the Newman- 
Penrose formalism. 



5.1. Initial data 



An important aspect in the solution in the time-domain of the perturbation equations Mil 
and ( I32> is the specification of the initial data, namely the determination of the Regge- 
Wheeler and Zerilli functions and "I^ ' at a given initial time and in the presence of the 
external sources of perturbations. The standard approach in this case is that of exploiting the 
1 + 1 splitting of spacetime and to solve the perturbed Hamiltonian and momentum-constraint 
equations as a set of coupled ordinary differential equations having the matter perturbations 
as sources. 

In the case of even-parity perturbations, in particular, the Hamiltonian constraint takes 
the form lElHol 



2 - 5M 



2M 



Or 



2M 



X 



1 



X 



-8ttT, 



on 



(45) 



and the momentum constraint is instead given by 



r - 2M 



dr*k 



2M 2 + rA(r - 2M) ■ ( 2M\ 4(r - 2M) - Ar 
-k + ( 1 ) — x 



2r 4 



4r 3 



r 



d r ,X = 



4tt 



2M 



4rd r „T*™ + - [AT? m + 4 (2r - 3M) T^] \ , (46) 



where k = dtk and \ = 9tX- 

Clearly, equations ( I45> and ( I46> are not sufficient to determine all of the unknowns: i.e. 
(k, k) and (x, X) an d m i s degree of indeterminacy reflects the arbitrary amount of gravitational 
radiation that can be present at the initial time with the same matter-sources. To alleviate at 
least in part this problem, the initial data is chosen so as to be best adapted to the physical 
problem under consideration. To minimize the initial gravitational-wave content, for instance, 
it is possible to set x = = x (conformally-flat approximation) and to use eqs. ( 145 \ and J46I 
to compute k and k. Alternatively, it is possible to consider time-symmetric initial data for 
which k = = x an d compute x an d k from the constraint equations [59 1. Finally, it is 
possible to specify a solution for x an d x an d use the constrain equations to compute k and 
fc[58|. 

Similarly, in the case of odd-parity perturbations, we can exploit the linearized Einstein 
equations [cf equations (18) and (20) of |39|] 



g CD V D [4r(V [A r)fc c] + 2r 2 V [c fc A] ] - (A - 2)fc y 



-167rL £ r 



(47) 



where 



L 



ivn 



V A k' 



16ttL 



Cm 



A(A - 1) J sin0 



tin 



sin 



2 n ^33 ) X-t m 



(48) 



(49) 
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Tedious but straightforward algebra allows to select two of the three equations d47b and J48l > 
as constraint equations for the four unknowns (fc , ko) and (fci, fci) 

d rM + ^T d rM ~ ^u k i ~ e 2b dr,k + ^-e 2b k = l6ird r ,L em (50) 
_a . 2(r-3M) Q , /2M A\ 1 , 2Af , 2 • 

d rM + ~ 2 d rt k + ^— --s -^k + - ^/C 

where Ua = c?t^A- As for the even-parity pertubations, the constraint equations J50i and ( 15 1> 
can be solved as a set of coupled ordinary differential equations to compute two of the four 
unknowns when the other two are specified according to the problem under consideration. 

As a final remark we recall that an equivalent and alternative approach to the specification 
of the initial data consist in setting = = and in solving the Regge-Wheeler 

and Zerilli equations as ordinary differential equations 1151 1161 . In the case of even-parity 
perturbations this is equivalent to setting fc = = x (cf equation (I44l ~) and thus considering 
time-symmetric initial data. 



5.2. Frequency domain and QNMs 

We conclude this Section by recalling that the investigation of the Regge-Wheeler and Zerilli 
equations (1211 — d32l> is often made in the frequency domain, where the the properties of 
Fourier transform can be exploited and the oscillation properties of the spacetime are easier 
to interpret. What is usually done in this case is to assume a harmonic time dependence in 
the perturbations, i.e. \I>( e ) ~ exp(«u; n i) and ~ exp(iui n t) where u n is the complex 
oscillation frequency of the n-th mode, so that the Regge-Wheeler and Zerilli equations (12 1> - 
d32t take respectively the form 

" + (u>l - V}° } ) *(°) = S<°> , (52) 

(V c >) " + (u* - V e ic) ) * (c) = , (53) 

where a 'prime' indicates a total derivative with respect to r* . Equations (1521 — d53l can be 
solved as an eigenvalue problem with boundary conditions such that the perturbations behave 
as pure outgoing-waves at spatial infinity 

$( c ) ~ cxp{-ioj n r*) for r* — > +oo . (54) 

and as pure ingoing-waves at the event horizon 

$r(°) ) $(°) ^ exp(zw n r») for r* -> -oo , (55) 

In the absence of source terms, the modes solutions of equations d52t and (1531 are referred to 
as the quasi-normal modes (QNMs) of the black hole, with the real part of uj n representing the 
oscillation frequency and with the imaginary part of lj^ 1 representing the damping time of 
the oscillation. A complete discussion of the properties and astrophysical relevance of these 
modes can be found in Q] |2 [5J and here we only recall the most important properties for a 
Schwarzschild black hole: 

• All of the QNMs have positive imaginary parts and represent therefore damped modes. 
Stated differently, a Schwarzschild black hole is linearly stable against perturbations. 
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• The damping time of QNMs depends linearly on the mass of the black hole and is shorter 
for higher-order modes. 

• The excitation of a black hole is referred to as its "ringing" and the late-time dynamics 
of these perturbations (i.e. the 'tail' of the ringing) can be described with a power-law 
representing the envelope of the various decaying QNMs. 

• The QNMs are isospectral, i.e. odd and even-parity perturbations have the same complex 
eigenfrequencies. This is due to the uniqueness in which a black hole can react to a 
perturbation and it is not true for a relativistic star. 

6. Asymptotic expressions and gravitational waves 

In the previous Section we have reviewed the derivation of the equations describing the 
evolution of perturbations of nonrotating black holes induced by a non-zero stress-energy 
tensor. These perturbations have been assumed to be generic in nature, needing to satisfy 
only the condition of having a mass-energy much smaller than that of the black hole. The 
solution of these equations with suitable initial conditions completely specifies the reaction 
of the black hole to the perturbations and this is essentially represented by the emission of 
gravitational waves. 

As mentioned in Section |2j the importance of the gauge-invariant variables used so far 
is that they are directly related to the amplitude and energy of the gravitational-wave signal 
measured at large distances. The purpose of this Section is to review the steps necessary to 
obtain the relations between the master functions for the odd and even-parity perturbations and 
the 'plus' and 'cross' polarization amplitudes h + , h x of a gravitational wave in a transverse 
and traceless (TT) gauge. In practice, and following the guidelines tracked in 8 1511161 . we will 
derive an expression for the perturbation metric h equivalent to that obtained in the standard 
TT-gauge on a Minkowski spacetime and relate it to the odd and even-parity master functions 
*(°)and*( e ). 

To obtain this result a number of conditions need to be met. Firstly, we need to evaluate 
each multipole of the decomposed metric perturbations in the tetrad e of stationary observers 
in the background Schwarzschild spacetime, i.e. hps, = e^e^h^, where e is diagonal with 
components = {e fc , e~ b , r _1 , (rsin^) -1 } and where the indices p refer to the locally 
'flat' coordinates. Secondly, all of the quantities need to be evaluated far away from the source 
(i.e. in the "wave-zone") and in the so-called radiation gauge. In practice, this amounts to 
requiring that components hgg, htt and hg^ are functions of the type fit — r)/r (i.e. they 
are outgoing spherical waves), while all the other components have a more rapid decay of 
C(l/r 2 ). Thirdly, we need to impose the condition that the metric is traceless modulo higher 
order terms, i.e. h§@ + = + 0(l/r 2 ). 

6.1. Asymptotic expressions from odd-parity perturbations 

We first consider odd-parity perturbations and recall that from the radiation-gauge conditions 
and since for large r the metric asymptotes that of a flat spacetime e b ~ e~ b ~ 1, we have 
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where the m indices have been omitted for clarity. Similarly, since h?f) — 2hr 2 VgSg ~ 
0(l/r) we deduce that h ~ 0(r) so that the only components of the metric having wave-like 
properties at large r are 

*-5(»»-»S)-?(^*-^)+ (p). ™ 

and where we note that since h has the dimensions of a length squared, both h + and h x are, 
as expected, dimensionless. 

Next, we need to relate the perturbation h to the odd-parity master function To do 
so, we follow the procedure outlined in 1 15 1, and note that (cf equation (111-20) |j) 

d t h=(l-^d r (r¥°))+h i ° ) , (60) 

Equation (I60> represents one of the Hamilton equations as derived by Moncrief in a 
Hamiltonian formulation of perturbation equations |14| (see |Appendix C| for details). We 
also note that the radiation-gauge conditions on h and imply that ^/(°) ~ O(l), i.e., in 
the wave-zone has the dimensions of a length, behaves as an outgoing-wave, but it does 
not depend explicitly on r. 

Exploiting now the outgoing-wave behaviour of h at large distances we can write 

d t h = -d r h + O (^) , (61) 

so that asymptotically equation i60\ simply becomes 

d r h = -d r (r*<°>) +0 (j) > < 62 > 
and its integration yields 



(63) 

As a result, the ' +' and ' x ' polarization amplitudes of the gravitational wave can be calculated 
from eqs. J59I — J58i as 

= (W^) + (?(l) , (64) 

(65) 

It is now interesting to note that expressions J64> and J65i can be written in a compact 
form using the s = — 2 spin-weighted spherical harmonics 



sin 



r sin \H 



so that expressions (I64> and (1651 can be combined into a single complex expression given by 



|| We recall that in the notation of 1151 i/> = A(A — 2)^(°\ and the multipoles in 1151 are related to ours as 
ha = ha = —2h, ho = h °^ and hi = 
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where, for clarity, we have explicitly restored the multipole indices t, m. 

An important clarification is now needed. We have discussed in the previous Section 
that the odd-parity metric perturbations are sometimes expressed in terms of the odd-parity 
Moncrief function Q(°> (cf equation J26H and it is not unusual to find in the literature the 
gravitational wave amplitudes expressed in terms of this quantity. However, great care must 
be paid to the asymptotic relation between the master function and the gravitational- 
wave amplitudes and, indeed, this is sometimes source of confusion |63 64 1. To clarify this 
point, we recall that the derivation of the asymptotic relation between and h proceeds as 
in a way similar to the one discussed above. In the radiation gauge and at large distances from 
the black hole, we can use relation J61t in the definition J26b with h,2 = —2h, so that 



)< o) ~ -d t h + o(- 



r \ r 



(68) 



which is also a dimensionless quantity. Since h ~ O(r), the function does not depend 
on r at leading order and equation d68l can be integrated to give 

/,(/) - ^ Q^(t')dt' + o(^j + const. , (69) 



r 

where the integration constant can be defined as 



const. = lim ~ 0(1) , (70) 



i- 



and it can be set to zero in the case of asymptotically flat metric perturbations (h = 0) at earlier 
times. Combining now expressions (1391 . J40i and ( I69K the gravitational-wave amplitudes in 
the two polarizations and with the new master function read 



Comparing expressions J67l > and (17 1 i we note that while 'J'^ -' and are both solutions 
of the Regge-Wheeler equation d2 U . they yield two different asymptotic expressions for the 
gravitational-wave amplitudes. This difference, which is consistent with equation ( I27> when 
evaluated in a an asymptotic region of the spacetime where L\ m = 0, is subtle but important 
and, as mentioned above, it has led to some inconsistencies in the literature both for the 
determination of the asymptotic gravitational-wave amplitudes and for the energy losses. This 
will be further discussed in sections |63]and|6!4| 



6.2. Asymptotic expressions from even-parity perturbations 

A calculation conceptually analogous to the one carried out in sections l6~T1 leads to the relation 
between the gravitational-wave amplitude and the even-parity master function. In particular, 
after projecting equation along the stationary tetrad, the asymptotic wave amplitudes in 
the two polarization states are 

& - \ OS - »S) - f ( w - *&r) - 1 • 
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so that we essentially need to relate the metric perturbation G with the even-parity function 
\j/( c ) Firstly, it is easy to realize that the even-parity metric projected onto the tetrad, h^l is 
such that 



H 2 ~O^J, and h^-O^-j, (74) 

so that the terms proportional to these multipoles are of higher order for large r and can be 
neglected. Furthermore, from the traceless condition 



4]+ h tl = + °{^) • (75) 
we obtain an asymptotic relation between K and G of the type 

V e V 9 Y em + U(2K-GA)F ta ~oU , (76) 

sin 6 J V r J 

where we have used the definition (|5j to derive the right-hand-side of expression (I76> . As 
a result, the asymptotic relation between the two components of the even-parity part of the 
perturbation metric simply reads 

K~^G + o(\) . (77) 



1 



Using now the definitions (I28t -(l29ll. we have that asymptotically 

ki ~ | G + rd r G + O (78) 

k 2 ~ - \ (K + rd r K) ~ — (G + rd r G) + O , (79) 

and their linear combination d29t becomes 

gi ~ ^ A(A-2)+offl . (80) 
Finally, the asymptotic gauge-invariant even-parity master function reads 

<gW~ a~~t / a ^ T ~-rG + ofiV (81) 

A [r(A - 2) + 6M] 2 W 

so that, modulo higher-order terms, the even-parity gravitational-wave amplitudes measured 
by a distant observer can be written in the compact form 



»S?-'»PL-;V^i»fi- I * , w*) +0 (?)- <82 > 

6.3. Asymptotic general expressions 

It is often convenient to combine the expressions for the asymptotic gravitational wave 
amplitudes due to odd and even-parity perturbations into a single expression of the type 
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or, equivalently 

(84) 

where we have defined h+ = + h+ and h x = h$ + h>x • Note that X m = for any 
value of I, so that in the case of axisymmetry the gravitational wave signal is proportional to 
l¥ ffl only. 

It is also useful to underline that while expression 1831 resembles the corresponding 
expression (10) of |63|, it is indeed very different. Firstly, because in [63 1 the Moncrief 
function is adopted for the odd-parity part of the perturbations and hence, modulo a 
normalisation factor, the function \&(°J appearing there corresponds to our function ( 
cf expression &26\ ). Secondly, because with this choice for the odd-parity perturbations a 
time derivative is needed in the asymptotic expression for the gravitational-wave amplitudes 
(cf the discussion in the derivation of equation (17 ID . As a result, expression (10) of |63| 
(which is also missing the distinction between the real and imaginary parts) should really be 
replaced by our expression (1841 . A similar use of the Moncrief function for the odd-parity 
part is present also in 1 651 1671 1581 . where it is employed to calculate the gravitational-wave 
content of numerically simulated spacetimes. 

6.4. Energy and angular momentum losses 

Using the expressions derived so far we can now estimate the energy and angular momentum 
losses due to gravitational waves propagating outwards to spatial infinity. More specifically, 
this can be done by using expression J83i and the definition of Isaacson's stress-energy 
pseudo-tensor t^ v for the gravitational-wave field h propagating in the curved background 
field g and in a Lorentz gauge 1691 1701 

1> = ^- (VnhapVuh *) , (85) 

where the brackets (. . .) refer to a spatial average over a length much larger than the typical 
gravitational wavelength |69, 11|. The averaged expression (JsHJ is gauge-invariant |69| 
and holds in the 'limit of high frequency' (or short-wave approximation), i.e. whenever 
the wavelength of the gravitational-wave field is small when compared to the local radius 
of curvature of the background geometry. In practice, gravitational radiation from isolated 
systems is of high frequency whenever it is far enough away from its source. 

Expression J85l > accounts for the amount of energy and momentum carried by the 
gravitational wave over a certain region of spacetime, but since we are interested in the 
energy flux as measured by an inertial observer, we need to project the pseudo-tensor on 
the observer's locally orthonormal tetrad, where it becomes 

r A£ > = ^ (dfihtfdfih**) , (86) 

with hjxo = hp^o — \hrjp,{>. As a result, the energy per unit time and angle carried by the 
gravitational waves and measured by a stationary observer at large distance from the black 
hole is given by 
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where the total derivative is made with respect to the asymptotic observer's time. Integrating 
( I87l i over the solid angle, the total power emitted in gravitational waves is then given by 



dE _ 1 (£_ - 

~dt ~ 167T f"^ JF- 



■2)1 



16tt 



2)! 



d^ 



(<0 



(if 



dV 



(e) 



f/f 



A 



A 



(88) 



(89) 



where expression ( I89l l was first presented in 115111 61 . 

For the same reasons discussed in the previous section, these expressions need to be 
suitably modified when the energy losses are expressed in terms of the odd-parity Moncrief 
function Qw, in which case they assume the form 

1 (£ + 2)1 ( d^ (c) 2 



dE 
~dt 



16tt 



(£-2)\ 



Cm 



dt 



O (o) 



(90) 



Besides having arbitrary angular dependence, the perturbations can also have non- 
diagonal terms and thus account for a non-zero net angular momentum, part of which can 
be lost to gravitational radiation. In this case, the angular momentum flux carried away in 
the form of gravitational waves can be calculated in terms of the r0 component of the stress 
energy tensor J86I . In particular, using spherical coordinates and assuming that the rotation is 
parametrised by the angle cj), we have 



d 2 J 



r 

167T 



drhfj^hgg + drh^d^h^ 



64> u 4>' 



(91) 



dtdVL 32tt 

Since the metric components with the radiation-gauge condition behave like outgoing 
spherical waves and since hgg = h + and hgi = h x , the angular momentum carried away 
in the form of gravitational waves ( 19 U is then expressed as 

■dthxd+ht) , (92) 



d 2 J r- 

= -— (dth+d^ 

l07T 



dtdD, 

Proceeding in a way similar to the one followed in the calculation of the emitted power, the 
total angular momentum lost per unit time to gravitational wave reads |49 1 



dJ 
~dt 



1 (1 + 2) 

Lm 



dV 



(<0 



dt 



(•2)- + *2(.SQ 



(£-2) 

or, using the Moncrief master function i26\ for the odd-parity perturbations (6^ 



dJ 
~dt 



> lm — 

16tt f-^ 



2)! 



(£-2)\ 



d^ 
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Q 



(°) 



(93) 



(94) 



To conclude, we report the expression for the energy spectrum dE/duo, which is readily 
calculated from equation (18 8 1 after performing the Fourier transform of the odd and even- 
parity master functions, i.e. 

dE _ 1 ^{£ + 2)1 
dF ~ 16^2 2^ (£-2)1 U 



(e) 



(95) 
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where we have indicated with f(u>, r) the Fourier transform of the timeseries f(t, r). As for 
the emitted power, also the energy spectrum will have a different expression if the odd-parity 
Moncrief function is used and in this case one obtains 

dE 1 v (l + 2)! / ~ (e) 2 (o) 2X 

We note that expression ( I96> for the energy spectrum agrees with the ones given in 1 53 56 1 for 
the cases of odd and even-parity perturbations, respectively and with the complete expression 
in 1 52 1 for both odd and even-parity multipoles. On the other hand, the energy spectrum 
reported in 1 63 1 suffers of the same inconsistencies discussed above for equation J84b . so that 
equation (27) in |63 1 should effectively be replaced by equation ( I96l l. 

We conclude this ection by noting that a new formalism has been recently proposed to 
compute the energy and angular momentum fluxes absorbed by a Schwarzschild (and Kerr) 
black hole in a form which is suitable for time-domain computations |62 49 1. Although the 
effect of the absorbed flux is generally small in astrophysical situations, it can be important in 
the secular dynamics in black- hole spacetimes, as discussed in |52|. Referring the interested 
reader to |62 49 1 for more details, here we simply recall that the calculation of the absorbed 
energy flux cannot be made with the present choice of Schwarzschild coordinates (in which 
the ingoing flux at the horizon is suppressed by the infinite redshift) and must therefore 
be performed using an advanced time-coordinate. Furthermore, Isaacson's prescription for 
the stress-energy pseudo-tensor J85i also ceases to be valid in the vicinity of the horizon, 
where the short-wave approximation breaks down and an alternative approach needs to be 
implemented 1621 1491 . 



6.5. A commonly used convention 

We conclude this section on asymptotic expressions by discussing a rather popular choice 
for the gauge-invariant master functions and that has found successful application in the 
extraction of the gravitational-wave content of numerically simulated spacetimes 1481 l34l 
l35l 1381 . Furthermore, the convention discussed below has been implemented in the 
Cactus computational toolkit 1361 1371 . a diffused and freely available infrastructure for the 
numerical solution of the Einstein equations 1711 . Numerous tests and applications of this 
implementation have been performed over the years and we refer the reader to I36ll37ll72ll73l 
for examples both in vacuum and non-vacuum spacetimes. 

The reference work for this convention in the one by Abrahams and Price (1996) |48|, 
although a similar approach for the even-parity part of the perturbations was also adopted in 
previous works 1441 1471 . We first note that the coefficients Co, c\ and C2 introduced in |48| 
are related simply to the multipolar coefficients of the odd-parity part introduced in section|3] 
More specifically, c 2 = — 2h — h 2 , c = and C\ — hi and it is then easy to realise 
that the odd and even-parity master functions Q^ m and Q\ m defined in |48] are related to the 
master functions discussed so far through the simple algebraic expressions 



2(1- 
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> (97) 



Qn,, = \i z rr J ^r*%> (98) 

so that the asymptotic expression for the gravitational-wave amplitudes in the two 
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polarizations are given by 
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(99) 



-y • (100) 



Also expressions ( I99> and J100I can be combined into a single one 



h+ - ih x = (Qtm - 1 / 

V2r f-f V J-, 



(t')dt' 



Imi 



O 



(101) 



which closely resembles expression J84b and that in its compactness highlights the advantage 
of the normalisation j97M98>. Also very compact is the expression for the emitted power 
that, with this convention, simply reads 

2 \ 



dE _ 1 ^ 



Lm 



dQ 



t-n i 



dt 



(102) 



We note that expression (I102> corrects equation (24) of |48|, where a time derivative of the 
odd-parity perturbations is present, leading to an obvious dimensional inconsistency. This 
typo was subsequently corrected in 1531 . 



7. Conclusions 



Black hole perturbation theory has been and will continue to be a powerful tool with 
continuous developments and multiple applications. In this paper we have briefly reviewed 
the theory of metric perturbations of Schwarzschild black hole, concentrating in particular on 
its gauge-invariant formulations. Special care has been paid to "filter" the technical details in 
the mathematical apparatus that may obscure the important results and to "steer" the reader 
through the numerous conventions and notations that have accompanied the development 
of the formalism over the years, pointing out the misprints and errors that can be found 
in the literature. Among the aims of this review is that of providing the reader with a set 
of expressions in the two most common formulations and that can be readily used for the 
calculation of the odd and even-parity perturbations of a Schwarzschild spacetime in the 
presence of generic matter-sources. In addition, we have reported the asymptotic expressions 
for the gravitational-wave amplitudes and losses driven by these generic perturbations. The 
latter can find useful application in the calculation of the gravitational-wave content of 
numerically simulated spacetimes. 

With suitable gauge choices, much of of what presented here for a Schwarzschild black 
hole could be extended to a Kerr background and this will be the subject of future work. 
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Appendix A. Some explicit expressions 



In the general expressions given in the main text the tensor harmonic appear in terms of 
covariant derivatives. While compact, such expressions need to be 'unwrapped' in practice 
and in what follows we provide their explicit expressions. 

VgV Y em = d}Y lm , (A.l) 

V V^y £m = dlY lm + sin 6 cos 8d e Y £m , (A.2) 
V V fl r £m = d^Y tm - cot 6d^Y lm , (A3) 
for the even-parity ones, and 

Si m = -^—d 4> Y lm , Si m = sm6d e Y em , (A.4) 

sin o v 

V 9 S e e m = - -i- {8^Y lm - cot 6d^Y lm ) , (A.5) 

V ^ m = sin 6 (8^Y £m - cot 6d^Y lm ) , (A.6) 
VgSf 1 = sin6d 2 Y lm , (A.7) 

WsSi 711 = - -^— n dlY tm - cos 6d 6 Y em . (A.8) 
sin 6 v 

for the vector and tensor odd-parity part. 

We also recall that the orthogonality relations between scalar, vector and tensor 
harmonics imply 



/ 
/ 
/ 



<Kl Y em Y t'm' = 5ee*S mm > , (A.9) 

^ 7 cd (d c y/ m )(d d w) = M u >s mm > , (a.io) 

dSlZ^Zti = A{A ~ 2) 6 u ,5 mm , , (All) 



while the orthogonality relations between the odd-parity axial vectors S c and tensor V^S^) 
follow straightforwardly from the above equations, giving 



/ 
/ 



dni cd s *em S ejm' = M u ,6 mm ,, (A.12) 
dnW (c S* d f n ^ c S d J m , = 2A(A - l)S u ,S mm , . (A. 13) 



Appendix B. A useful example: point-particles 

In the spirit of providing a useful example of a matter-source generating odd and even- 
parity perturbations of a Schwarzschild spacetime, we now a particular problem that has been 
carefully studied over the years: that of a point-particle orbiting around a black hole. Since 
this problem has been solved in many different cases and using very different approaches, it 
can serve as a useful test for codes handling more complex matter-sources. Because of this, 
we report below the relevant components of the odd and even-parity decomposition of the 
particle's energy-momentum tensor. 

For a particle of mass fi -C M, specific energy E and angular momentum L, the 
equations of motion in the equatorial plane (6 — tt/2) for the coordinates R(t) and $(£) 
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are 
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where R = dR/dt, $ = d^/dt and where the energy-momentum tensor reads 
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As a result, the odd-parity components of the stress energy tensor are^f 
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while the even-parity components assume the form 
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Appendix C. Variational principle for odd-parity perturbations 

In Sect. |4] we have discussed that it is possible to choose two different representations for 
the odd-parity perturbations, namely expressions J20> for the CPM convention and expression 

% Note that the spherical harmonics and their derivatives are all evaluated at <f> = <J>(t) and 9 = w/2. 
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fl26l > for the RW convention. An intimate connection links the two master functions and this 
can be appreciated by using the variational formalism employed by Moncrief in 1141 . To 
highlight this, let us recall that two perturbation functions k\ , k 2 can be introduced for the 
odd-parity perturbations 



,(°) 



d r h 2 



-2h. 



(CI) 



where k\ is gauge-invariant (cf equation ( I18» but k 2 is not. While this may seem just 
another possible choice for the odd-parity perturbations, it points out that it is now possible to 
introduce the quantities tti and tt 2 as conjugate momenta of k\ and k 2 in the Hamiltonian 

2" 



drH = \ 
A 



dr 



-A 



1 

r 

dr 



2r(r - 2M) 
(A -2) 



7T 2 



A-2 



1 



1 n 1 

-O r 7Ti TTi 

2 r 



drh^ir 2 



(C2) 



In addition to the constraint equation ir 2 = (which must be enforced at initial time and is 
then conserved since k 2 is cyclic), the other Hamilton equations are 



SHn 



dth 



SHrr 



(C3) 



<57Ti ' " " 5tt 2 ' 8k\ 

where S( )/S( ) refers to a functional differentiation. The first and third of eqs. 1C.3I have 
explicit expressions 
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A(A - 2) 
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2M 
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5 r 2 (r 2 7n) 
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, (C4) 



(C5) 



and the Regge-Wheeler equation in the RWM convention can be obtained by time- 
differentiating equation dC.4> . replacing dt~Ki through equation dC.5> and introducing the 
quantity 

2Af 



1 - 



(C6) 



to obtain 

a 2 o(°) - a 2 ow + f 1 - — 

Equation dC.7> coincides with equation (4.20) of 1 14 1 with Q = Q^°\ 

Remarkably, the Regge-Wheeler equation in the CPM convention can be obtained also 
through the Hamilton equations in terms of the conjugate moment to ki, i.e. rw\. More 
specifically, it is sufficient to differentiate in time equation dC.5> and to replace dtki in 
equation (IC.4> to obtain 

2M\ / A 6M\ 
1 — - -=- )rw 1 = , (C.8) 



Equation 1C.7I coincides with equation (IT 15) of 1 15 1 with 



t 

A 



d t h { ° ] - r 2 d r 



(C9) 
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Note also that the second of Hamilton equations ( IC.3I can be used to derive 

2 /, 2M\ ,.,..> > „, (o) 



dth 2 



1 



-) d r (r 2 7n) - 2^° 



(CIO) 



A(A-2) 

which coincides with equation (111-20) of Cunningham et al. 1151 and has been used to obtain 
the asymptotic relation between $ and the gravitational wave amplitude (cf equation 

Finally, we underline that the relation \21\ between the two master functions and 
Qv°) can be obtained in coordinate independent form also in the formalism introduced by 
Martin-Garcia and Gundlach |39|. More specifically, exploiting the odd-parity perturbations 
equations J47i as 



g CD V D [4r(V [A r)fc c] + 2r 2 V [c k A] ] n A -(A- 2)k A n A = -16nL e ^n A (C.ll) 

and the gauge-invariant and coordinate independent definitions of and (i.e. 
equations dl9l and ( I26> of the main text), straightforward calculations then yield 

Q( o) 16tt 



i>(°) - 



where * (o) 



u A V A ^ {o) and r' 



A-2 rn La 



(C12) 



n Vat are 'frame derivatives' in terms of the non- 



coordinate basis vectors [u , n ) of M , such that — uau = hati — 1, 9ab = 
—uaub + tiatib and €ab = tlaub — ua^b- In the case of Schwarzschild coordinates, 
n A = (0, e~ b ), u A = (e b , 0) and equation 1C. 121 coincides with equation ( I27> . 

Appendix D. Angular pattern functions 

Finally, this section is devoted to a list of the expressions of the angular functions W lm and 
X im for the first values of t and to. While these expressions can be easily reproduced with 
straightforward algebra, they are tedious to derive and hard to find in the literature. 
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X 32 = ie i2 V sin 9 cos 20 . 
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q / c 

W 40 = —= (5 + 7 cos 26») sin 2 , W 41 = e^-W - (7cos 20 - 1) sin 20 , (D.7) 

8-y/7T 8 V 71" 
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